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MULTIPLE-MODE NONLINEAR FREE AND FORCED VIBRATIONS 
OF BEAMS USING FINITE ELEMENT METHOD 

'if 'ic'if 

Chuh Mei and Kamolphan Decha-Umphai 
Old Dominion University, Norfolk, Virginia 

ABSTRACT 

Multiple-mode nonlinear free and forced of beam has been analyzed 
by the finite element method. The geometric nonlinearity is 
investigated. Inplane (longitudinal) displacement and inertia { ID I ) are 
also considered in the formulation. Harmonic force matrix is derived 
and explained. Nonlinear free vibration can be simply treated as a 

special case of the general forced vibration by setting the harmonic 
force matrix equal to zero. The effect of the higher modes is more 

pronounced for the clamped supported beam than the simply supported 
one. Beam without inplane displacement and inertia (IDI) yields more 
effect of the higher modes than the one with inplane displacement and 
inertia. The effects of inplane displacement and inertia are to reduce 
nonlinearity. [For beams with end supports restrained from axial 
movement (immovable case), only the hardening type nonlinearity is 
observed]. However, beams of small slenderness ratio (L/R = 20) with 
movable end supports, the softening type nonlinearity is found. The 
concentrated force case yields a more severe response than the uniformly 
distributed force case. Finite element results are in good agreements 
with the solutions of simple elliptic response, harmonic balance method, 
Runge-Kutta method and experiment. 

* Professor, Department of Mechanical Engineering and Mechanics 
** Graduate Student 
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Chapter 1 
INTRODUCTION 

In modern engineering, with its continuous refinement of instrumen- 
tation, its improved computational capabilities, and the high precision 
tolerances, the theory of nonlinear vibrations is gaining more and more 
practical meaning. Although it is known that linear vibrations provide 
no more than a first order approximation of an actual situation, they 
are sufficient for some practical and engineering purposes. The linear 
theory is inadequate, however, if the vibration of an elastic body 
involves amplitudes that are not very small, as assumed in the linear 
theory. Nonlinear vibration approach leads to completely new phenomena 
which are not possible in linear systems. For example, the dependence 
of frequency, or period of vibration, on amplitude cannot in principle 
be handled by using linear analyses. In such cases, nonlinear theory 
must be used to obtain more accurate results and to explain new 
phenomena. The steadily increasing demand for more realistic models of 
structural responses has resulted in solution techniques to deal with 
nonlinear structural problems. Apart from some few exceptions, it is 
generally not possible to provide analytical closed-form solutions for 
the differential equations occurring in nonlinear vibrations of 
structures. Naturally, a numerical solution may be obtained when the 
motion corresponding to certain initial and boundary conditions is to be 
determined. Because of the advance in electronics, modern digital 
computers have been of great value in solving nonlinear problems. 

In general, nonlinearities in structural mechanics problems can 
arise in several ways. When material behavior is nonlinear, the 
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generalized Hooke's law is no longer valid. This type of nonlinearity 
is called "material" or "physical" nonlinearity. Alternatively, 
material behavior can be assumed to be linear, but structural 
deformation can become large and cause nonlinear strain-displacement 
relations. Deformation of a structural member can also reach a 
magnitude that does not overstrain the material; in such a case, 
curvature of the deformed median line can no longer be expressed by a 
linear equation. Problems involving large structural deformation are 
called "geometrically" nonlinear problems. A combination of material 
nonlinearity and geometric nonlinearity is also possible. 

Any large-amplitude deflection of a beam which is restrained 
axially at its two ends, results in some midplane stretching. This 
stretching must be accounted for by using nonlinear strain-displacement 
relationship (geometric nonlinearity). The nonlinear equation of motion 


describing this 

situation 

had 

been the 

basis of 

a number of 

investigations. 

Most of 

these 

works are 

based on 

a single-mode 


approach. Woinowsky-Krieger 1 considered the effect of an axial force on 
the vibration of hinged bars. The vibration of an extensible bar, 
carrying no transverse load and having the ends fixed at the supports, 
caused an axial tensile force with a period equal to the half-period of 
the transverse vibration of the bar. The elliptic function was used to 
produce the relation of frequency and amplitude of vibration. Eringen^ 
studied the nonlinear free vibrations of elastic bars having immovable 
hinged ends. The solution was accomplished by the use of the 
perturbation method. The ratio of nonlinear period over linear period 

and axial stress were shown against initial deflection. Only the 
hardening type nonlinearity was found. Burgreen^ studied the nonlinear 
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free vibrations of a pin-ended column whose ends were pinned to points 
fixed in space. This imposed the condition of a constant end distance 
instead of the usual theoretical assumption that the axial load in the 
column remained constant along the beam length. The elliptic function 
was performed to get the exact solution. He also found that the 
frequency was dependent upon the amplitude of vibration, the effect of 
the amplitude of vibration becoming more pronounced as the Euler load 
was approached in which the classical linear theory yielded the 
frequency of vibration as zero. Woodall^ considered the nonlinear free 
vibrations of a thin elastic beam. In his formulation, a fixed inertial 
reference frame and a Lagrangian description of the motion were 
employed. By assuming the motion to be inextensional and, at the same 
time, admitting the existence of a resultant normal force acting on each 
cross-section of the beam, a system of governing equations was 
obtained. The solutions of the simply-supported beam were obtained by 
using three methods: the finite difference method, perturbation 

technique and Galerkin weighted residual method. For the particular 
example considered in his paper, the finite difference solution appeared 
to be stable, even for oscillations involving angular rotations at the 
ends of the beam of the order of magnitude of 80°. Furthermore, he 
found that the Galerkin approximate solution was in closer agreement 
with the finite difference solution than the perturbation solution. 
Raju et al. 5 studied free flexural vibrations of a simply-supported beam 
when a compatible longitudinal or inplane mode was coupled with the 
fundamental flexural mode. The Rayleigh-Ritz method was employed. The 
results showed that the effects of longitudinal displacement and inertia 
were to reduce the nonlinearity in the flexural frequency-ampl i tude 
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relationship. Tseng and Dugundji 6 investigated a straight beam with 
fixed ends excited by a periodic motion at its supporting base in a 
direction normal to the beam span. By using Galerkin's method, the 
governing partial differential equation was reduced to the well-known 
Duffing equation. The harmonic balance method was applied to solve the 
Duffing equation. Pandalai 7 investigated the case of straight beams, 
irrespective of the boundary conditions, the nonlinear was found to be 
of the hardening type. He further concluded that only the hardening 
type can be existed. Later, Atluri 8 showed that there were some cases 
for which the softening type nonlinearity is possible. He investigated 
the large amplitude transverse vibration of a hinged bar with one end of 
the beam free to move longitudinally. The equation was solved by the 
perturbation procedure of multiple-time scales. The calculated results 
showed that the predominant nonlinear effect was due to longitudinal 
inertia which was of the softening Duffing type. This result was in 
contrast to the earlier analyses where a hardening nonlinearity had been 
predicted when the only nonlinearity considered was the effect of 
average midplane stretching due to the out-of-plane deflection. 

Through the foregoing studies, the general features of the 
nonlinear response of beams under harmonic excitation seem to have been 
largely clarified. However, most of the actual calculations are based 
on a one-term approximation for the spatial function while the 
interactions between the modes with amplitude are not addressed. 
Corresponding to this phenomena, the effect of multiple modes on beam 
response is needed. Because of the complexity in the formulation, there 
are very few investigations on multiple-mode analysis in the 
litera ture. 
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McDonald^ apparently was the first to consider modal 

interactions. The considered investigation was the vibration of a 
uniform beam with hinged ends which were restrained. The beam was 

subjected to a concentrated lateral force at the mid-span and then 

released from rest at the deflected position. The nonlinear effect in 

this investigation was produced by the axial stretching of the beam. By 

assuming a multiple-mode expansion corresponding to the deflected 
position, the elliptic function procedure was performed to evaluate the 
coefficients related to the participation of each mode. Bennett and 
Eisley 1 ® investigated the steady-state free and forced responses and 
stability for large amplitude nonlinear vibration of a beam with clamped 
ends. The general equations for the response and stability were 
derived. By applying Galerkin's method, a set of nonlinear ordinary 

differential equations was obtained. The solution of forced vibration 
was evaluated by the method of harmonic balance. Later, Bennett 11 
considered the problem involved the ul traharmonic response of a simply 
supported beam. Tseng and Dugundji 1 ^ also used a multiple-mode 
expansion in considering the forced response of a clamped beam about its 
buckled configuration. The buckled beam was excited by the harmonic 
motion at its supporting base. By using Galerkin's method, the 
governing partial differential equation reduced to a modified Duffing 
equation which was solved by the harmonic balance method. Srinivasan 1 ^ 
solved free and forced responses of beams undergoing moderately large 
amplitude steady-state oscillations by the averaging method of Ritz. 
The application of Ritz's method to solve the governing nonlinear 
partial differential equation, yielded nonlinear algebraic equations 
instead of nonlinear ordinary differential equations. To solve these 
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nonlinear algebraic equations, the Newton's method or bigradient matrix 
method had to be employed. The method was shown by assuming the first 
two symmetric modes of the linear system for the deflection of beam. It 
was clear that the method yielded as many simultaneous nonlinear 
algebraic equations as the number of modes included. Nayfeh et al. 1 ^ 
proposed a numerical-perturbation method for the determination of 
nonlinear forced response of beams. The deflection curve of the beam 
was represented with a multiple-mode expansion in terms of the linear 
modes. Then the temporal problem was solved by the method of multiple 
scales, and internal resonances was also considered. Van Dooren and 
Bouc 15 considered the nonlinear transverse vibrations of a uniform beam 
with ends restrained and forced transversely by a two-mode function 
which was harmonic in time. A simply supported beam was considered by 
the two-mode approach. Approximate solutions were found by using 
Urabe's numerical method applied to Galerkin's procedure and by an 
analytical harmonic balance-perturbation method. The existence of sub- 
harmonic response of order 1/3 and harmonic response in the sub-harmonic 
region of the forcing function was proved. Takahashi 1 ® analyzed the 
inextensible clamped-free and free-free beams by using Galerkin's method 
and the harmonic balance method. Yamaki and Mori 17 investigated 
nonlinear forced vibrations of a clamped beam under uniformly 
distributed periodic lateral loading with the effects of both initial 
deflection and initial axial displacement taken into consideration. The 
problem was first reduced to that of a finite degree-of-freedom system 
with the Galerkin procedure, the steady-state solutions which were 
obtained by applying the harmonic balance method. Actual calculation 
was carried out for the three degree-of-freedom system with symmetric 
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modes. Yamaki et al. also performed experiments to compare to the 
analytic results. The test results were reported in the root-mean- 
square of deflection instead of the actual deflection of the beam. 

Generally, the classical approach to solve nonlinear vibrations of 
a beam is to start with the so-called assumed mode shape. By employing 
the Galerkin's method, the governing nonlinear partial differential 
equation of motion is reduced to a system of nonlinear ordinary 
differential equations. The elliptic function, perturbation method, or 

numerical methods, can then be employed to solve the problem. 

In practice, many optimum or minimum-weight designed structures are 
complex. Because of the versatility of the finite element method, it is 
more suitable to use this method to analyze complex structures. Mei 1 ^ 
investigated nonlinear vibration of beams by matrix displacement 
method. Nonlinear free vibrations of various boundary conditions were 
investigated and good agreements were obtained between the finite 
element method and other numerical methods. Rao et al.^ studied the 
large amplitude free oscillations of beams and orthotropic circular 
plates. Their finite element formulation was based on an appropriate 
linearization of the nonlinear strain-displacement relations. Simply 
supported and clamped beams were investigated. Comparison of their 
results with the earlier work confirmed the reliability and 
effectiveness of the linearization of the strain-displacement 
relations. Reddy and Singh^ A investigated the large-deflection analysis 
of thin elastic curved beams by conventional and mixed finite element 
methods. The conventional finite element method was based on the total 
potential energy expression, whereas the mixed method was based on a 
Reissner-type variational statement and involved the bending moments and 
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deflections as primary dependent variables. From their result, it 
appeared that in general the mixed method yielded more accurate 
results. Recently, Mei and Decha-Umphai 22 * 23,24 developed the harmonic 
force matrix for solving nonlinear forced vibrations of beams and plates 

or 

by finite element method. By employing the linearizing method* 13 and so- 
called "single-mode" approach, the frequency-amplitude-force relations 
were obtained. There are very few attempts to solve nonlinear 
vibrations of beam using finite element method with multiple-mode 
expansion. Busby and Weingarten 26 investigated beams of simply 
supported and clamped boundary conditions. Their finite element 
technique was performed only to obtain the nonlinear differential 
equations of the straight beam and the method of averaging was then used 
to obtain an approximate solution. Cheung and Lau 27 investigated the 
two-mode nonlinear vibration of beams. The essence of their method 
could be regarded as an incremental harmonic balance method associated 
with a finite strip procedure in the time-space domain. Unfortunately, 
the works reported in References 26 and 27 could not be exactly 
classified as finite element method. 

The accuracy of the theoretical predictions would not be completed 
unless the experimental studies had been compared. Bennett and Eisley 10 
performed the experiment of a clamped beam to compare with their 
theoretical results. Tseng and Dugundji^ conducted the experiment of a 
straight beam with fixed ends which was excited by the periodic motion 
of its supporting base in a direction normal to the beam span. All of 
these two experiments were found to compare favorably with the corre- 
sponding theoretical predictions. However, these experiments were 
carried out with special kinds of excitations, eg. supporting base 
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excitation. Yamaki et al. 18 performed the experimental studies of a 
clamped beam under a uniformly distributed periodic load. Besides the 
reasonable agreement with the theoretical predictions, 17 their 
experimental results seem to provide effective data facilitating further 
theoretical analyses. 

It is clear that a substantial amount of literature exists on 
nonlinear vibrations of beams. Sathyamoorthy^ 8 *^ published two 
excellent survey articles, one which dealt with literature concerning 
classical nonlinear methods, and another with the finite element 
method. 

Through the foregoing studies, the general features of the 
nonlinear free and forced vibrations of beams seem to have been largely 
clarified. However, modern structures are complex and the more accurate 
theoretical predictions are preferred, therefore, the multiple-mode 
approach has to be considered in the formulation of nonlinear vibration 
problems. Since the evolution of digital computers, the finite element 
method has become widely used to solve many types of complex 
structures. It is the purpose of this research to extend the finite 

element method to multiple-mode nonlinear vibrations. 

In this thesis research, multiple-mode nonlinear free and forced 
vibrations of beams using finite element method are presented. Both 
out-of-plane deflection and inplane displacement are included in the 
formulation. The classical method is provided in Chap. 2 to review the 
concepts of nonlinear vibrations. In Chap. 3, the finite element 
formulation is presented in detail. The nonlinear stiffness matrices 

pc 

using the linearizing method^ 3 are derived. The harmonic force matrix 
for multiple-mode approach is derived for both uniformly distributed and 
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concentrated forces. In Chap. 4, the solution procedure of multiple- 
mode nonlinear vibrations by finite element method is presented. In 
Chap. 5, beams with various out-of-plane and inplane boundary conditions 
are investigated. The definition of inplane boundary condition is also 
explained. The relations of frequency and amplitude for various 
boundary conditions and loads are tabulated and plotted. Results are 
also compared to other classical or experimental results whenever 
available. The convergence criteria are also studied. Chapter 6 gives 
the concluding remarks. 
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Chapter 2 

CLASSICAL FORMULATIONS 

In this chapter, the mathematical formulations of the classical 
method are expressed. Equations of motion for an isotropic beam are 
presented. The characteristic equations are provided by performing the 
Galerkin's method. For single-mode approach, the frequency-amplitude 
relation is shown in a simple closed-form solution which could not be 
done for the multiple-mode approach. For multiple-mode approach, the 
general formulation is shown to provide for better understanding in the 
multiple-mode formulation. All of these formulations are performed in 
detail . 


2.1 Classical Method 

In general, nonlinear forced vibration of a beam is solved by the 
assumed mode method. By employing the Galerkin's method, the assumed 
mode shape is substituted into the governing equation of motion, then 
the integration is performed over spatial domain. The characteristic 
equations are obtained. The numerical time integration or other 
numerical method is needed for solving the frequency-amplitude 
relations. Both single and multiple-mode approaches are shown in the 
formulations. Because of the simplicity of single-mode approach, this 
approach is shown first to provide some concepts which later is extended 
to the multiple-mode formulation. 
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2.2 Single-Mode Approach 

The classical method for single-mode approach is straight 
forward. The frequency-amplitude relation exists in a simple closed- 
form relation which could not be done for multiple-mode approach. The 
formulation for single-mode approach is as follows. 

Assume that a uniform beam with cross sectional area A, moment of 
inertia I and length L is subjected to a uniformly distributed periodic 
lateral load F(x,t) as shown in Fig. 1. In this figure, a clamped beam 
is shown for simplicity. The deflection is denoted by w(x,t). With the 
assumption that axial displacements at both ends are zero (immovable 
inplane displacement), the basic governing equation of motion for the 
nonlinear bending forced vibration of a beam is found to be 

pA w + El w, xxxx - N w, xx - F(x,t) = 0 (2.1) 

where 

EA ^ 2 

N«^J <w, x ) dx (2.2) 

o 

where p and E are mass density and Young's modulus, respectively. The 
subscript x following a comma denotes differentiation with respect to x 
and a dot denotes differentiation with respect to time t. 

For single-mode approach, the deflection is assumed as 

w(x,t) = 4> ( x ) g(t) (2.3) 

where <|>(x) represents the normal mode of linear free vibration which 
satisfies the boundary conditions for the case in consideration, and the 
modal amplitude g(t) is an unknown function of time. 
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Using the foregoing expression for w, the Galerkin's method can be 
applied to Eq. (2.1), which leads to 
L 

/ [Eq. (2.1)] 4> ( x ) dx = Residual = 0. (2.4) 

o 

Equation (2.4) can be expressed as 

_ L 2 .. L 

[pA / <t> dx] g + [El / 4>, xxxx 4> dx] g. 


L L 

- [N / <t>, <t> dx] g - F(t) / (|» dx * 0. (2.5) 


Equation (2.5) can be rewritten in the form as 


pA g + [El / i>, xxxx $ dx] g // <|> 2 dx 

o o 


L L o L L , 

- [N / 4>, xx 4> dx] g/J 4> dx - F(t) / dx// <J> dx = 0. (2.6) 

O O 0 0 

By substituting Eq. (2.3) into Eq. (2.2), the inplane force N can be 
expressed as 

N ■ [ f£ J U, x ) 2 dx] g 2 . (2.7) 

0 

By using this value of N, Eq. (2.7), in conjunction with Eq. (2.6), the 
characteristic equation is obtained in the form as 

m g + kg + k g 3 - c F(t) = 0 (2.8) 


where m, k and k are the mass, linear stiffness and nonlinear stiffness 
terms, respectively. The values of m, k and k can be expressed as 
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m * p A (2.9) 

L L ? 

k * El / <t>» xxxx <J> dx/J 4> dx (2.10) 

o o 

EA 2 ^ *■ 2 

k = J (♦, ) dx / ♦, <t> dx/ / <D dx. (2.11) 

cL X XX 

0 0 0 

The value of constant c is expressed as 

L L 9 

c = / <t> dx// <D dx. (2.12) 

0 0 

The characteristic equation, Eq. (2.8), can be written as 

m g + kg + k - P(t) = 0 (2.13) 

where P(t) is the force term which is expressed as 

P(t) - c F(t) (2.14) 

The characteristic equation, Eq. (2.13), can be solved for the 
frequency-amplitude-force relation in a closed-form solution by various 
approximate methods. It should be pointed out that many of these 
approximate methods yield the same closed-form solution which can be 
classified as a standard form. Because of the simplicity of the 
following method, it is chosen to be shown in depth for the derivation 
of the closed-form solution. This method is based on omitting the 

higher harmonic term. 

To obtain the frequency-amplitude-force relation in closed-form 
solution, g(t) is specified as 

g(t) = A cos (wt) (2.15) 

where A is the amplitude of deflection of the beam and u> is nonlinear 
frequency. The nonlinear term, g 3 (t), is expressed as 
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g 3 (t) = A 3 cos 3 (cot). (2.16) 

By expanding cos 3 (cot) , it can be written as 

3 3 1 

cos (cot) = - cos (cot) + -r cos (3 cot). (2.17) 

4 4 

By neglecting the higher harmonic term, cos (3 cot) , the approximation 
of cos'* (cot) can be written as 

cos'* (cot) = 7 cos (cot). (2.18) 

4 


By substituting Eq. (2.18) into Eq. (2.16), the nonlinear term is 
expressed as, 

g 3 ( t ) =2 A 3 cos (cot). (2.19) 

4 


The force term, P(t), is written as, 

P(t) = P Q cos (cot) 

and by Eq. (2.14), 


and 


C F 


o 


F ( t) = F q cos (cot) 


( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 


where F Q is force intensity which has the dimension as force per unit 
length. 

By substituting Eqs. (2.15), (2.19) and (2.21) into Eq. (2.13), the 
frequency-amplitude-force relation can be shown as 

2 3-3 

- co m A cos (cot) + k A cos(cot) + — k A cos(cot) 


- P Q cos (cot) = 0. 


(2.23) 


Equation (2.23) can be rewritten in the closed-form solution as 
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w 2 = a> L 2 + | 0 A 2 - P o /A (2.24) 

where 

u> L 2 = k/m (2.25) 

p » k/m . (2.26) 

Equation (2.24) is the standard closed-form relation between frequency 
and amplitude and force. w L is defined as linear frequency which has 
the dimension as radian per second. 

2.3 Multiple-Mode Approach 

The classical method for multiple-mode approach follows the same 
path as single-mode approach. But the multiple-mode approach could not 
yield the simple closed-form frequency-amplitude-force relation, thus, 
numerical integration or other approximate method is needed. 

The general formulation of multiple-mode approach is as follows, 

n 

w(x,t) * 2 <t>.(x) g. (t) (2.27) 

i=l 1 1 

where ^(x) is the i-th normalized linear mode shape with maximum of 
unity, and g^t) is the time function. The total number of modes in 
consideration is n. 

The time function g.j(t) can be expressed as 

g^t) = H. cos(wt) (2.28) 

where is the amplitude of i-th mode shape. By using the assumed 
expression of w, Eq. (2.27), the Galerkin's method can be applied to Eq. 
(2.1) as follows: 
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L 

/ [Eq. (2.1)] <j>. (x) dx * Residual = 0 (2.29) 

o 1 

i 1)2. « • i n. 

By substituting Eq. (2.27) into Eq. (2.2), the axial force N is 
expressed as 

N = -|r / [2 <(>. g.]^ dx. (2.30) 

o i=l 1,x 1 


By substituting Eq. (2.30) into Eq. (2.29), the following set of 
nonlinear ordinary differential equations are obtained: 


m i 


k i 9i 


n n n 
+ I E Z 
j=l r=l s=l 


k. . g.q q 

Tjrs a j a r a s 


where 


= Pj(t) i = 1,2, •••, n 


m. = p A 


k f = EI I *1,xxxx »f dx/ -f 4 dx 


/ <!> . <t> dx / $ <{». 

ri « J» x r » x S,XX 1 

i, _ tA 0 0 

ijrs 2L C ~ 


dx 


/ 4> i dx 
o 


P- (t) = / F(x,t) <|». dx// <b. 2 dx. 

3 A ' ^ » 


(2.31) 

(2.32) 

(2.33) 


(2.34) 


(2.35) 


The set of nonlinear ordinary differential equations shown by Eq. (2.31) 
is the set of characteristic equations for nonlinear forced vibration. 
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In general, the number of characteristic equations is the same as the 
number of mode shapes which are included in the assumed deflection 
shape, Eq. (2.27). These characteristic equations are highly nonlinear 
and coupled, thus, it is very tedious to solve for the steady-state 
solution. One way to obtain the frequency-amplitude-force relation is 
by performing numerical integration, e.g., Runge-Kutta method. 

To clarify the method, a two-mode approach is performed as follows: 

2 

w(x,t) = r (x) g - ( t) . (2.36) 

1=1 

The deflection shape w(x,t) is assumed as the combination of two 
modes. In this case, n is equal two. The value of N, Eq. (2.30), can be 
expressed as 


N =# [(/ *1% M 


9 ? + 


(/ 4 y dx) 


4 


+ 12 I *l.x *2,x dx) 


(2.37) 


By using Eqs. (2.36) and (2.37), the set of characteristic equations, 
Eq. (2.31), can be obtained as 


m 


1 9 1 + k l 9 1 + k Ull 9 1 + ^ 1112 9 1 9 2 + k 1 122 9 1 9 2 


+ k 1222 g 2 3 - Pj_ (t) 


(2.38) 


m 2 g 2 + k 2 g 2 + *2111 9 1 3 + *2112 g l 2 g 2 + *2122 g l g 2 2 


+ k 2222 9 2 = p 2 (t) 


(2.39) 
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where 


ir> 1 = P A 


k l ■ EI ! ♦l.xxxx *1 ix, l *1 2 dx 


EA , L 2 L 


k im = 2L 1 *l,x dX J *l,xx *1 dx// *1 dx 
0 0 0 


r !» I ! , , l ,, L 2 

k 1112 - 2L 1 *l,x dX 1 »2.xx *1 dX// * 1 dx 
0 0 0 


EA L 


*1122 S 7C f Q *l,x *2.x dx ^ *2, xx *1 dx/ ^ Q *1 dx 


EA r L 2 L 


k 1222 = 2T J *2,x dX 1 *2,xx *1 dx// *1 dx 
0 0 0 


P^(t) * / F(x,t) 4>j dx// $ dx 
0 0 




4> 


2,xxxx 


L 


<t> 2 dx// 

o 



dx 


21 


EA L 2 L 2 

k = — f < b dx / d> <b dx/f d> dx 

2111 2L J v l,x J l,xx *2 J ? 2 
oo o 


- EA f L 2 J f L J /f L 2 ^ 

k = — / d> dx j d) <b dx/ J <t> dx 
2112 2L J 9 l,x J 9 2,xx ? 2 J 9 2 
0 0 0 


- EA ^ ** *■ 2 

k 2122 = 2L / V *2,x ** ! *2,xx‘2 dx/J *2 dX 
0 0 0 


EA L 2 L *■ 2 

k * — f d> dx / d> d> dx/ / d> dx. 
2222 2L J V 2,x J V 2,xx 9 2 /J V 2 
oo o 


(2.40) 


The steady-state response for a two-mode nonlinear vibration 
problem can be obtained by employing the Runge-Kutta method for solving 
the nonlinear coupled characteristic equations, Eqs. (2.38) and (2.39). 
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Chapter 3 

FINITE ELEMENT FORMULATION 

In this chapter, the formulation of a finite element method is 
presented. The expressions for strain-displacement relation, kinetic 
energy and potential energy are provided. The linearizing function for 
deriving nonlinear stiffness matrices is defined. The inplane 

displacement and lateral deflection are expanded in the cubic order. 
Furthermore, the derivation of harmonic force matrix is performed. 

3.1 Strain and Curvature-Displacement Relations 

The strain-displacement relation of a beam element shown in Fig. 2 


is given by 



. 1 2 

e * u. + - w. 

(3.1) 

X X Z X 

where e is the strain in x-di recti on of beam, u 

A 

is the inplane 

displacement and w 

is the lateral deflection. 

The curvature- 

displacement relation 

is defined as 



^x ” w ’xx 

(3.2) 


where <|» x is the curvature. The total strain, c , is defined as 

12 /.ox 

£ = u, + - w, - z w, . (3.3) 

x 2 x xx 




24 


3.2 Kinetic and Strain Energies 

The kinetic energy of a beam element is given by 

T = jj p J (u 1 2 + w 2 ) dx (3.4) 

o 

where p is mass per unit length and l is element length. The kinetic 
energy can be separated into two parts as 

T = T $ + T b . (3.5) 

The kinetic energy due to mid-plane stretching T s is defined as 

l 

T $ ■ g P / u dx (3.6) 

0 

and the kinetic energy due to bending T b is defined as 

1 

T b ■ g p / w dx * (3.7) 

o 

The strain energy of a beam element is given by 

. I 

U - 5 J (E A t + E I 4 T ) dx (3.8) 

w XX 

0 

where E is Young's modulus, A is the cross-sectional area of beam, and I 
is the moment of inertia. By substituting Eqs. (3.1) and (3.2) into Eq. 
(3.8), the relation between strain energy and displacements can be 
written as 

l 2 

U = | J E 'a [u.* + 2 U , X (J «.*> - (| -.*) ] dx 
0 

1 A 2 
+ - / E I w, dx. 

2 xx 

o 


(3.9) 


The strain energy U can be separated into two parts, namely, linear 
strain energy U L and nonlinear strain energy U NL . 

U = U L + U N| _. (3.10) 

The linear strain energy for Eq. (3.10) is expressed as 

JL X 

U L = \ I EAu,[d x + 2 J El w,^ x dx. (3.11) 

0 0 

Similarly, the nonlinear strain energy for Eq. (3.10) is expressed as 

JL 2 

U ML * | / E A [2 u, x <i w, x ) * (i w, x ) ] dx. (3.12) 

0 


3.3 Displacement Functions 


The displacement functions are chosen to be in cubic order as 
follows: 


and 


w = a^ + a 2 x + a 3 x + a 4 x~ 


u = a 5 + a g x + a 7 x + a g x' 


(3.13) 


(3.14) 


The generalized coordinates aj, , a g can be written in vector form 

as 


{a} = [a^ a 2 a g a 4 a g a g a 7 a g ] 


(3.15) 


and the element nodal displacement. Fig. 2, are defined as 
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where 0. and a., are the differentiation with respect to x of w^ and u n - 
at node i, respectively. The generalized coordinates in Eq. (3.15) can 
be written as 

(a) = [T] (6) (3.17) 

where [T] is a transformation matrix. 

The displacements u and w including their derivatives can be 
expressed in terms of element nodal displacements as 


u = [0 0 0 0 1 x x 2 x 3 ] (a) 

= [B] [T] {6} (3.18) 

u. = [0 0 0 0 0 1 2x 3x 2 ] {a} 

A 

= [C] [T] {6} (3.19) 

w * [1 x x 2 x 3 0 0 0 0] (a) 

= [D] [T] {6} (3.20) 

w, =[01 2x 3x 2 0 0 0 0] {a} 

A 

= [E] [T] {6} (3.21) 

w, xx = [0 0 2 6x 0 0 0 0] (a) 


= [G] [T] {6} 


(3.22) 
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3.4 Linearizing Function 


The linearizing function is defined as 5,25 

f = 4 w, 

2 ’x 


(3.23) 


By substituting Eq. (3.21) into (3.23), the linearizing function can be 
expressed as 

f CE] CT] (6>. (3.24) 


This linearizing function f is assumed to be constant in each element. 
The main advantage of using linearizing function is shown in the 
following sections. 


3.5 Element Equations of Motion for Nonlinear 
Free Vibration 

The nonlinear free vibration is the backbone of the investigation 
of steady-state response. The equations of motion can be derived by 
applying Lagrange's equations as follows. 

From Eqs. (3.6) and (3.18), the expression of kinetic energy due to 
mid-plane stretching T s can be written in matrix form as 

T = i (6} T [T] T / p [B] T [B] dx [T] {6} 

S c 

0 

= 4 {*> T &" 3 «> (3.25) 

2 s 

where [m $ ] is the element mass matrix due to inplane displacement 
function and is expressed as 
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T X Y 

Cm s 3 = [T] / P CB3 CB] dx [T]. (3.26) 

Similarly, from Eqs. (3.7) and (3.20) the expression of kinetic energy 
due to bending can be written as 

T « 4 <$} T [m. ] {6} (3.27) 

b i b 

where Cm b ] is the element mass matrix due to bending and is expressed as 

[m b ] = [T] T / p [D] T [D] dx CT]. (3.28) 


The total kinetic energy T, Eq. (3.5), can be expressed in matrix form 
as 

T = “ (6} T [m] {5} (3.29) 

where [m] is the element mass matrix and is expressed as 


Cm] = Cm $ ] + Cm b ] 


(3.30) 


The linear strain energy U L , Eq. (3.11), can be expressed in matrix 
form as 

U L * l - {6} T Ck L ] (6> (3.31) 

where Ck b ] is the element linear stiffness matrix. The element linear 
stiffness matrix Ck|_] can be separated into u-part Ck LS ] and w-part 
Ck Lb ] as 

Ck L ] = Ck L$ ] + Ck Lb ] (3.32) 

where 

Ck L$ ] = CT] T / E A CC] T CC] dx CT] 


(3.33) 
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T A T 

[k L() ] = [IT / El CG] 1 [G] dx [T]. (3.34) 

o 

The nonlinear strain energy U NL , Eq. (3.12), can be written by 
using the linearizing function f, Eq. (3.23), as 

1 X 

U.„ = i / E A [2f u, w, + f w, ] dx (3.35) 

NL C XX X 

0 

By substituting Eqs. (3.19) and (3.21) into Eq. (3.35), the nonlinear 
strain energy can be expressed in matrix form as 

v ■ \ {6)T c V {5> <3,36) 

where [k NL ] is the element nonlinear stiffness matrix and is expressed 
as 

X 

Ck NL ] = [T] T / {E A f [CC] T [E] + CE] T [C3] 

0 

+ E A f 2 [E] T [E]} dx [T]. (3.37) 

The linearizing function f which is a constant for each element 
transforms the element nonlinear stiffness matrix [k^] into the linear 
form as in Eq. (3.37). The strain energy U, Eq. (3.10), can be written 
as 

U = | {6} T [k] {6} (3.38) 

where [k] is the element stiffness matrix and is expressed as 


[k] = [k L ] + Ck NL 3. 


(3.39) 



The application of Lagrange's equation leads to the element nonlinear 
free vibration equation of motion as 
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Cm] {6} + [[k^] + [k^]] (6) = 0. (3.40) 

3.6 Element Harmonic Force Matrix 

In Chap. 2, the equation of motion for nonlinear forced vibration 
in classical approach is expressed in Eq. (2.13) as 

m g + kg + k g 3 = P(t) (3.41) 

when the forcing function P(t) is a simple harmonic P Q cos (u>t), an 

approximate solution of Eq. (3.41) is in the form of Eq. (2.24) as 

“ 2 = 2 + J M 2 - P /A. (3.42) 

L 4 0 

When the forcing function P(t) is a simple elliptic function and 
expressed as 

P(t) = B* A cn (\t, rj ) 

* B*g (3.43) 

where g is the Jacobian elliptic function and expressed as 

g - A cn (Xt.n) (3.44) 

★ 

where B is forcing amplitude factor, \ is circular frequency of the 
Jacobian elliptic function and n is modulus of the Jacobian elliptic 
function. 

By expanding the elliptic forcing function into Fourier series then 
comparing the order of the magnitude of the various components, Hsu 3 ® 
concluded that the harmonic forcing function P Q cos (wt) is the first 
order approximation of the elliptic forcing function B A cn (\t, n). 


Also, the first order appoximation of the elliptic response of Eq. 
(3.41) yields the same frequency-amplitude-force relation of Eq. (3.42) 
which is the perturbation solution (standard form, Eq. (2.24)). To 
obtain the exact elliptic response of Eq. (3.41), the forcing function 
P(t) in Eq. (3.43) is treated as a linear spring in the Duffing 
equation, Eq. (3.41). 

•• - 3 * 

mg + kg+kg = B g 


m g + (k - B*) g + k g 3 = 0. (3.45) 

★ it O 

This linear spring force B g possesses a potential energy as B g /2. 
Similarly, the potential energy of a beam due to the uniform harmonic 
forcing function F Q cos wt can thus be approximated by 


V = E 3 / w 2 dx (3.46) 

c o 

where the summation denotes the sum of all elements. 

To find the value of B*, the conclusion by Hsu 30 is needed (Eq. 
(3.43)) as follows 

P(t) = B* g 

= B A cn (xt,n) 


so. 


P„ COS(a)t) 
0 


~ * 


B A cn 


(3.47) 


By using Hsu's conclusion, the harmonic function cos (wt) is the first 
order approximation of the elliptic function cn (Xt.n), thus. 
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By using the relation of Eq. (3.48), substitutes into Eq. (3.47), the 
relation of B* can be expressed as 

B* = P q /A. (3.49) 

From Eq. (2.21), the relation of P Q is expressed as follows 

P o a c F o (3 - 50) 

where the constant c is expressed in Eq. (2.12). The constant c can be 

interpreted as the ratio of area under mode shape and square of mode 

shape. And the nonlinear mode shape in multiple-mode approach is 

assumed to be sum of linear modes, thus, the deflection w(x,t) can be 

written as * * 

w (x,t) = <t> (x) g (t) (3.51) 

where 4 (x) is normalized nonlinear mode and g*(t) is time function, 
thus, 

* 1 n 

4> (x) * - E A. 4 >. (x) (3.52) 


where Aj and <(> ^ ( x ) are the amplitude of i-th mode shape and the i-th 
normalized linear mode shape, respectively. By using Eq. (3.52) and the 
definition of c in Eq. (2.12), in general, the expression of c for 
simple harmonic force is expressed as 


(Load elements) 



L *2 


/ ( 4> ) dx 


o 


c 


dx 


(3.53) 
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Thus, for the uniformly distributed force over the entire beam, the 
value of c is expressed as 


r * 

) <f> dx 

c = -j? • 

/ U*) 2 dx 

o 

In the finite element method, the concentrated force can be represented 
by the distributed force over a small length, thus, the value c for 
concentrated can be expressed as 

b * 

/ dx 

c = (3.55) 

/ U*) 2 dx 

o 

where a and b are coordinates of beam under that distributed force. 

By using the Eq. (3.46), the element harmonic force matrix can be 
derived from 

* X 

V = •£— / W 2 dx (3.56) 

o 

where B* is expressed in Eq. (3.49). By substituting Eq. (3.20) into 
(3.56), the potential energy due to the uniform harmonic force can be 
expressed as 

* * 

V = | (6} T [T] T / [0] T [D] dx [T] {6} 
w 0 

= \ (6} T [h] {6} (3.57) 

where [h] is the element harmonic force matrix and expressed as 
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[h] = B* [T] T / [D] T [D] dx [T]. (3.58) 

o 

By comparing Eq. (3.58) to Eq. (3.28) and using the expression of B* and 

P Q , Eqs. (3.49) and (3.50), the harmonic force matrix can be written as 

c F 

[h] = — — [m, ] (3.59) 

A p b 

where c is a constant expressed in Eqs. (3.53), (3.54) or (3.55). It 
should be noted that the harmonic force matrix is a symmetrical matrix. 

3.7 Element Equations of Motion for Nonlinear 
Forced Vibration 

The application of Lagrange's equation leads to the equations of 
motion for nonlinear forced vibration of a beam element under harmonic 
forcing function. By using Eqs. (3.29), (3.38) and (3.57), the 

equations of motion are expressed as 

Cm] {5} + [[k^] + [k NL ] - [h]] {6} = 0 (3.60) 

where [m] is the element mass matrix in Eq. (3.30), [k L ] is the element 
linear stiffness matrix in Eq. (3.32), [k^] is the element nonlinear 
stiffness matrix in Eq. (3.37), and [h] is the element harmonic force 
matrix in Eq. (3.59). It should be noted that all these matrices are 
symme trie. 
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Chapter 4 

SOLUTIONS PROCEDURES 

In this chapter, the solution procedures of finite element method 
are explained in detail. The solution procedures are divided into two 
major parts, namely, the small deflection part which is the linear 
solution, and the large amplitude part which is the nonlinear 
solution. Each part consists of minor steps which evaluate the element 
and system matrices. The iterative process for the large amplitude part 
is also explained. The convergence criteria are also provided and 
convergence characteristics are shown in figures. Finally, the computer 
flow-chart of the solution procedures is provided at the end of this 
chapter. 


4.1 Small Deflection Solution 
(Linear Solution) 

Linear solution is performed to evaluate the linear eigenvalues 
(linear frequencies) and linear eigenvectors (linear mode shapes). The 
combination of these linear mode shapes and amplitude for each mode is 
represented for multiple-mode approach which is later used in the 
nonlinear solution. 

The application of Lagrange's equation leads to the equation of 
motion for linear free vibration of a beam element as 

[m]{6 } + [k^]{5> = 0 


( 4 . 1 ) 
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where [m] and [k^] are the element mass and linear stiffness matrices as 
defined in Eqs. (3.30) and (3.32), respectively. The element mass 
matrices, [m $ ] and [m b ] as defined in Eq. (3.30), can be evaluated from 
Eqs. (3.26) and (3.28), respectively. After all element mass and 
stiffness matrices are known, the equations of motion for the system can 
be assembled as 

CM] (A) + [K l ] {A} = 0 (4.2) 

where [M] and [K L ] are the system mass and linear stiffness matrices, 
respectively. (A) is the eigenvector which consists of the vectors in 
term of w and w, , and u and u, . This eigenvector (A) can be expressed 

a A 

as 



where { 4>> is the eigenvector related to w and w >x> and {£} is the eigen- 
vector related to u and u, . By using the bending part of Eq. (4.2), 

A 

the system equations of motion for bending only can be expressed as 

[M b ] U) + CK Lb 3 {$} = 0 (4.4) 

where [M b ] and [K Lb ]are the bending system mass and linear stiffness 
matrices which are related to Eqs. (3.28) and (3.34), respectively. 
Equation (4.4) can be rewritten as 

“L CM b ] { \ } = CK Lb ] <*L> (4.5) 

where co^ is the linear frequency (eigenvalue) and {<j> b } is the corre- 
sponding linear mode shape normalized to unity at the middle of beam. 
It should be pointed out that for each linear frequency, there is a 
corresponding linear mode shape, for example, the first linear 
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frequency u> L1 corresponds to {<|> L1 }, o) L2 corresponds to and 

u Ln corres P° nds t0 

Through the foregoing procedure, the linear solution for free 
vibration of beam has been clarified. The linear frequencies; 

^Ll* u l 2* ***» ^Ln ancl the corresponding linear mode shapes; 
{4> l1 >, ^Ln^ have 136611 evaluate<1 * These linear mode shapes 

will be used later for multiple-mode approach in the next section. 

4.2 Large Amplitude Solution 
(Nonlinear Solution) 

In this section, the multiple-mode nonlinear free and forced 
vibrations of beam by finite element method are presented. The 

iterative procedure and the convergence criteria are also included. 

After the linear mode shapes have been determined (see Section 4.1 
Small Deflection Solution), the nonlinear solution is performed as 
follows. 

By using the first linear mode shape and the total amplitude at the 
middle of the beam, the deflection shape for the first iteration is 

expressed as 

{<(»} = A {<{>qK (4.6) 

By using Eq. (4.6), the deflection shape of each element {6} can be 

found. By using {6} and Eq. (3.24), the linearizing function ,f, for 
each element can be determined. After the linearizing function has been 
determined, the element nonlinear stiffness matrix [k kM ] for each 
element can be evaluated by using Eq. (3.37). The integration in 
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Eq. (3.37) can be performed by the extended Simpson's rule for 20 
interval s 32 . 

The element stiffness matrix [k] and element mass matrix [m] can be 
determined by using Eqs. (3.39) and (3.30), respectively. 

The element harmonic force matrix [h] can be determined by using 
Eq. (3.59). The constant c in Eq. (3.59) is expressed in Eqs. (3.53), 
(3.54) and (3.55). 

The application of Lagrange's equation leads to the element 
nonlinear forced vibration of beam under harmonic forcing function. 
This equation of motion is expressed as 

Cm] {6} + [[k] - Ch3 ] (6} = 0. (4.7) 

Equation (4.7) can be rewritten as 

Cm] {6} + [[k L ] + [k NL ] - [h]] {6} = 0. (4.8) 

After all the element mass, stiffness and harmonic force matrices 
have been determined for each element, the equations of motion for the 
system can be assembled as 

CM] (A) + [[K l ] + [K nl ] - [HD ] (A) = 0 (4.9) 

where CK^^] and [H] are the system nonlinear stiffness and system 
harmonic force matrices, respectively. The system equations of motion, 
Eq. (4.9), can be reduced into the terms of w amd w, x by using Guyan's 
reduction technique 31 and the definition of (A) is defined in Eq. 
(4.3). This reduced system equation of motion can be expressed as 


[RM] {(j)} + [RK] {<)>} = 0 


(4.10) 
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where [RM] is the reduced system mass matrix and [RK] is the reduced 
system stiffness matrix. It should be noted that [RK] consists of not 
only the reduced system linear and nonlinear stiffness matrices but also 
the reduced system harmonic force matrix. 

In the process of transformation Eq. (4.9) into Eq. (4.10) there 
exists the relation between the system deflection shape {A} and the 
reduced system deflection shape as 

{A} = [TRF] {*} (4.11) 


where [TRF] is the transformation matrix. 

By the definition of multiple-mode approach, the reduced system 
deflection shape { 4>} is the combination of the linear mode shapes 
{<j>^> in Eq. (4.5) and their amplitudes. This reduced system deflection 
shape {<t>} can be expressed as 


{*} a Aj {4> L1 > + A 2 {$ L2 } + •••• + A n {<t> Ln } 


(4.12) 


where A^ is the amplitude of the i-th mode, < - > is the linear mode 
shape of the i-th mode and n is the total number of modes. Equation 
(4.12) can be rewritten in a matrix form as 


or 


where 


f A ' 
A 1 


{$} = [(<t> L1 > U L2 > •••• 4| n }] 


4} = [<D Q ] (A 0 > 


Ln < . 


(4.13) 


(4.14) 


[<t> 0 ] = [{4> l1 } ^4>j_2 ^ ^Ln^ 
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<v = 



(4.15) 


By using Eq. (4.14), Eq. (4.10) can be transformed to the normal 
coordinates, {A q >, and is expressed as 


[<t> 0 ] T CRM] [<t Q ] {A q } + l> 0 ] T CRK] C <t> Q ] (A q > = 0 
or 

[ R M 1 { A } + [RK n ] {A > = 0. 

0 0 0 0 


(4.16) 

(4.17) 


Equation (4.17) is in the form of an eigenvalue problem which can be 
expressed as 

wjj L [RM q ] {A q } = [RK q ] { A q } (4.18) 

where is the nonlinear frequency. By solving the eigenvalue 

problem, Eq. (4.18), the nonlinear frequency and the eigenvector 
{A } can be determined. After the eigenvector {A } has been known, the 

U 1 U 

Ai 

amplitude ratios •»—; i*l,2,*»», n can be determined. 

A 1 

By the definition of multiple-mode approach, Eq. (4.12), the 
maximum deflection (amplitude) is the summation of amplitudes of all the 
modes. This can be expressed as 


or 


n 

A = 2 A 
i=l i 


A 

A n i 

— = 1+2 — 

A i=2 A 

1 1 


(4.19) 


(4.20) 


or 
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A = 
1 


n i 

l + z — 

i=2 A 


A i 1 

The ratio -r— is from Eq. (4.18) as mentioned before. 

A 1 

amplitude for i-th mode, A.,-, can be solved by 


(4.21) 


The value of 


A i 

A i " 


V 


i = 2,3, 


n. 


(4.22) 


Through the foregoing procedure, the first iteration has been 
completed. The next iteration starts by using 

{$} * ? A. U (4.23) 

1*1 1 Li 

By using Eq. (4.23) instead of Eq. (4.6), the next iteration is 
performed from that point onward to obtain the improved nonlinear 
frequency u> NL1 and i-th amplitude A.,-. This iterative process can now be 
repeated until a convergence criterion (Section 4.3) is satisfied. It 
should be noted that the nonlinear stiffness matrix Ck Nl _] and the 
harmonic matrix [h] are updated in each iterative process because of the 
changing of the values of A^. The flow-chart of this solution procedure 
is shown in Fig. 3. 


4.3 Convergence Criteria 

Three displacement convergence criteria proposed by Bergan and 
Clough"^ and a frequency convergence criterion are employed in the 
present study. The three displacement norms (criteria) are the modified 
absolute norm, the modified Euclidean norm and the maximum norm. The 
definitions of these three norms are given in the Appendix. The 
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Fig. 3 Computer Flow-Chart (Solution Procedures) 
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frequency norm is defined | ( Aco^^ ) j ) j I where (Au)^)j is the change 
in nonlinear frequency during the j-th iterative cycle. Fig. 4 shows a 
typical plot of these four norms versus number of iterations for a 
three-mode simply supported beam (L/R = 1010) without inplane 

displacement and inertia ( ID I ) subjected to a uniformly distributed 
harmonic force Fo = 0.001 N/rnm at A/R = 2.0. Fig. 5 shows a plot of 

these four norms versus number of iterations for a three-mode clamped 

beam (L/R = 1010) with inplane displacement and inertia (IDI) when both 
ends are immovable (u = 0 at x = 0 and L) and the beam is subjected 
to a uniformly distributed harmonic force Fo = 0.002 N/mm at 

A/R = 4.0. For all of the following results presented in this study, 
convergence is considered to be achieved whenever any one of the norms 
reaches a value of 10" 5 . 
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Chapter 5 

RESULTS AND DISCUSSIONS 

The fundamental frequency ratio (w = of multiple-mode 

nonlinear free and forced vibrations at various amplitudes for simply 
supported and clamped beams are reported in this chapter. Both 
immovable and movable inplane edges conditions are considered. Finite 
element results with and without inplane displacement and inertia (IDI) 
are given. The meaning of "without inplane dispalcement and inertia" is 
to neglect {£}, Eq. (4.3), completely from the formulation. The 
harmonic balance solution 17 , Runge-Kutta solution and experimental 
result 18 are also given for comparison with the finite element 
results. Because of symmetry, only a half of beam divided into twenty 
elements of equal length is considered herein. 

5.1 Boundary Conditions and Material Properties 

The transverse deflection boundary condition for simply supported 
beam is defined by letting the deflection (w) equal to zero at that 

location. For clamped supported beam, the transverse deflection 
boundary conditions are the deflection (w) and its associated 

slope (w, ) equaled to zero at that location. The inplane boundary 

A 

conditions are divided into two catagories; namely, immovable and 
movable inplane edge. The definition of immovable inplane edge is 

defined as the inplane displacement (u) at the boundary equaled to 
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zero. For the movable Inplane edge, the inplane displacement (u) is set 
free at the boundary. 

All the results presented here are based on the following material 
properties: 


mass density 

thickness 

width 

Young's modulus 
radius of gyration 


p = 26.6832 E-10 
h = 0.514 mm 
B = 26.0 mm 
E = 6.98 E+04 
R = 0.148379 mm 


N - sec 2 /mm 4 


N/mm 


2 


For nonlinear vibration, the effect of slenderness ratio (L/R) has 
some influence in the solutions. In this report, there are many 
slenderness ratios in use as shown in Table 1. These slenderness ratios 
are calculated by changing the beam length and keeping the cross- 
sectional area constant. Most of the results reported herein, is based 
on the 150 mm. beam length (L/R = 1010) otherwise specified. This 150 
mm. long beam (L/R = 1010) is the same dimension as the beam which 
Yamaki et al.*® performed the experiment. 


5.2 Improved Nonlinear Free Vibration 

The fundamental frequency ratios of ^ ree Vlbratl0n at 

various amplitudes (A/R) without inplane displacement and inertia (IDI) 
for clamped and simply supported beams (L/R = 1010) are shown in Table 2 
for the cases of single, two and three-mode method, respectively. The 
amplitude ratio for these beams are also provided in Table 3. Table 2 
shows that the more numbers of modes are used in the analysis, the less 
the frequency ratios will be, e.g., at A/R = 5.0, the three-mode 
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Table 1 Relations Between Slenderness Ratio 
(L/R) and Beam Length 


Slenderness 

Beam 

Ratio 

Length 

L/R 

L (mm.) 

1010 

150.000 

100 

14.840 

50 

7.420 

20 

2.968 
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Table 2 Frequency Ratios for Nonlinear Free Vibration of 
Clamped and Simply Supported Beams (L/R = 1010) 
without Inplane Displacement and Inertia ( I D I ) 


A/R 


Frequency Ratio, u> ,/u. . 

NL1 LI 



1 mode 

2 modes 

3 modes 


Finite 

Element 

Elliptic 

Solution 

Fini te 
Element 

Fini te 
Element 

Clamped Beam 

1.0 

1.0218 

1.0222 

1.0218 (2) a 

1.0217 (3) 

2.0 

1.0845 

1.0857 

1.0844 (3) 

1.0831 (4) 

3.0 

1.1817 

1.1831 

1.1814 (5) 

1.1757 (4) 

4.0 

1.3056 

1.3064 

1.3051 (6) 

1.2900 (6) 

5.0 

1.4495 

1.4488 

1.4490 (8) 

1.4188 (9) 

Simply Supported Beam 

1.0 

1.0897 

1.0892 

1.0888 (3) 

1.0888 (3) 

2.0 

1.3229 

1.3178 

1.3120 (5) 

1.3119 (5) 

3.0 

1.6394 

1.6257 

1.6030 (7) 

1.6022 (7) 

4.0 

2.0000 

1.9760 

1.9248 (11) 

1.9218 (11) 

5.0 

2.3848 

2.3501 

2.2624 (17) 

2.2549 (18) 

a. 

Number in brackets 
converged solution 

denotes the 
10" 5 

number of iterations 

to get a 
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Table 3 Amplitude Ratios for Nonlinear Free Vibration of Clamped 
and Simply Supported Beams (L/R = 1010) without Inplane 
Displacement and Inertia (IDI) 


A/R 


Amplitude Ratio 


2 modes 

3 modes 


Aj/A2 

Ai/Ag 

V a 3 

Clamped Beam 

1.0 

- 1054 

- 1056 

2164 

2.0 

- 280 

- 282 

559 

3.0 

- 137 

- 139 

261 

4.0 

- 87 

- 89 

157 

5.0 

- 64 

- 66 

109 


Simply Supported Beam 


1.0 

446 

446 

* 

2.0 

125 

125 

7652 

3.0 

66 

66 

1972 

4.0 

45 

44 

848 

5.0 

35 

34 

480 

* Number 

4 

is larger than 10 
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solution yields smaller frequency ratios than the two-mode solution. 
Table 3 shows the influence of the amplitudes of the higher modes, 
especially the amplitude of second mode. Table 4 shows the comparison 
of a two-mode response between the finite element method and the Runge- 
Kutta method for clamped and simply supported beams (L/R = 1010) without 
inplane displacement and inertia { ID I ) . This clearly demonstrates the 
remarkable agreement between finite element and Runge-Kutta solutions. 
These results are also shown in Figures 6 and 7 for the clamped and 
simply supported cases, respectively. 

The clamped and simply supported free vibration results 
(L/R = 1010) with inplane displacement and inertia (IDI) for both ends 
restrained from longitudinal movement (immovable case) are shown in 
Table 5. The amplitude ratios for this clamped beam are also shown in 
Table 6. The frequency-amplitude relationships for three-mode responses 
of these clamped and simply supported beams are also shown in Figure 
8. This figure clearly shows that the simply supported beam yields the 
higher nonlinearity than the clamped beam. 

The responses for both ends free to move longitudinally (movable 
case) are shown in Tables 7 and 8. Table 7 shows the frequency- 

amplitude relations for three-mode solution of clamped and simply 
supported beams with inplane deformation and inertia for the slenderness 
ratio of 100, 50 and 20. The resul t for the clamped case is also shown 
in Figure 9. From this figure, the high slenderness ratio beam 

(L/R = 100) yields almost none of nonlinearity. Conversely, the less 
slenderness ratio case (L/R = 20) leads to a situation that eventually 
exhibits slightly softening type nonlinearity. The softened type exists 
when the nonlinear frequency is less the linear frequency (a><1.0). 
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Table 4 Comparison Between Runge-Kutta Method and Finite 

Element Method for Two-Mode Nonlinear Free Vibration 
of Clamped and Simply Supported Beams (L/R = 1010) 
without Inplane Displacement and Inertia (IDI). 


A/R 

Frequency Ratio, 

“NLl^Ll 

% Difference* 


Finite Element 

Runge-Kutta 


Clamped Beam 

1.0 

1.0218 

1.0222 

0.04 

2.0 

1.0845 

1.0852 

0.07 

3.0 

1.1817 

1.1810 

0.03 

4.0 

1.3056 

1.3009 

0.32 

5.0 

1.4495 

1.4373 

0.81 


Simply Supported Beam 


1.0 

1.0888 

1.0888 

0.00 

2.0 

1.3120 

1.3135 

0.11 

3.0 

1.6030 

1.6115 

0.53 

4.0 

1.9248 

1.9467 

1.13 

5.0 

2.2624 

2.3023 

1.73 


* % Difference = 


(Runge-Kutta) - (Finite Element) „ inn 
CRunge-Kutta) u 




Fig. 6. Clamped beam W/0 IDI, 2 modes, L/R = 1010, free vibration finite element method vs 
Runge-Kutta method. 
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FREQUENCY RATIO <u Nu /u> u 

Free vibration, simply supported beam W/0 IDI, 2 inodes, L/R = 1010, finite element method vs, 
Runge-Kutta method. 
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Table 5 Frequency Ratios for Nonlinear Free Vibration of 
Clamped and Simply Supported Immovable Beams 
(L/R = 1010) with Inplane Displacement and Inertia 
( IDI ) 


Frequency Ratio, 


A/R 

1 mode 


2 modes 

3 modes 

Fini te 
Element 

Rayleigh 

Ritz 5 

Finite 

Element 

Fini te 
Element 

Clamped Beam 

1.0 

1.0149 

- 

1.0149 (2) a 

1.0149 (2) 

2.0 

1.0582 

- 

1.0581 (3) 

1.0581 (3) 

3.0 

1.1268 

- 

1.1264 (4) 

1.1259 (4) 

4.0 

1.2164 

- 

1.2151 (5) 

1.2140 (2) 

5.0 

1.3226 

• 

1.3202 (6) 

1.3176 (6) 


Simply Supported Beam 


1.0 

1.0607 

1.0607 b 

1.0607 (2) 

1.0607 (2) 

2.0 

1.2247 

1.2246 

1.2247 (2) 

1.2247 (2) 

3.0 

1.4577 

1.4573 

1.4577 (2) 

1.4577 (2) 

4.0 

1.7320 

1.7309 

1.7320 (2) 

1.7320 (2) 

5.0 

2.0310 

2.0289 

2.0310 (2) 

2.0310 (2) 


a. Number in brackets denotes the number of iterations to get a 
converged solution 10" 5 


b. L/R = 100 


Table 6 Amplitude Ratios for Nonlinear Free Vibration of 
Clamped Immovable Beams (L/R * 1010) with Inplane 
Displacement and Inertia (IDI) 


A/R 


Amplitude Ratio 


2 modes 

3 modes 


Aj/A2 

a x /a 2 

V*3 

1.0 

- 1354 

- 1354 

★ 

2.0 

- 346 

- 346 

2690 

3.0 

- 159 

- 159 

1202 

4.0 

- 94 

- 95 

678 

5.0 

- 64 

- 64 

442 


* 4 

Number is larger than 10 
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Table 7 Frequency Ratios for Three " Mo<le Nonlinear 

Free Vibration of Clamped and Simply Supported Movable 
Beams with Inplane Displacement and Inertia ( ID I ) for 
Different Slenderness Ratio (L/R) 


A/R 


Slenderness Ratio, L/R 


100 

50 

20 

Clamped Beam 

1.0 

.9999 (2) a 

.9997 (2) 

.9983 (2) 

2.0 

.9997 (2) 

.9989 (2) 

.9933 (2) 

3.0 

.9994 (2) 

.9976 (2) 

.9850 (3) 

4.0 

.9989 (2) 

.9957 (2) 

.9737 (3) 

5.0 

.9983 (2) 

.9933 (2) 

.9596 (3) 


Simply Supported Beam 


1.0 1.0000 (2) 

.9999 (2) 

.9993 (2) 

2.0 .9999 (2) 

.9996 (2) 

.9973 (2) 

3.0 .9998 (2) 

.9990 (2) 

.9938 (2) 

4.0 .9996 (2) 

.9982 (2) 

.9891 (2) 

5.0 .9993 (2) 

.9973 (2) 

.9832 (2) 

a. Number in brackets denotes the 

r- 

number of iterations 

to get a 


converged solution 10” 5 . 


Table 8 Amplitude Ratios for Three-Mode Nonlinear Free 
Vibration of Clamped and Simply Supported Move- 
able Beams (L/R = 20) with Inplane Displacement 
and Inertia (IDI) 


A/R 

Amplitude Ratio 


a^/a 2 

V*3 

Clamped beam 

1.0 

6410 

★ 

2.0 

1613 

* 

3.0 

724 

- 8914 

4.0 

414 

- 5152 

5.0 

270 

- 3413 


Simply Supported Beam 


1.0 

★ 

★ 

2.0 

★ 

★ 

3.0 

5178 

★ 

4.0 

2936 

* 

5.0 

1898 

★ 

* Number 

4 

is larger than 10 
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Atluri 8 also obtained the similar softening type in his investigation. 
For the simply supported case as shown in Table 7, it exhibits less 
influence of softening type than the clamped case. Table 8 shows the 
amplitude ratios for three-mode solution of clamped and simply supported 
movable beams (L/R = 20) with IDI. The comparison of the amplitude 
ratio between the movable clamped case (Table 8) and the immovable 
clamped case (Table 6) shows that the higher modes have more influence 
on the immovable case than the movable case. 

5.3 Nonlinear Response to Distributed 
Harmonic Force 

The responses of clamped and simply supported beams (L/R = 1010) 
without IDI are shown in Tables 9 and 10. Table 9 shows the frequency 

ratios for the cases of single, two and three-mode approaches. It 

should be noted that as the amplitude is increased, the more iteration 

is needed. The amplitude ratios for these beams is shown in Table 10. 

The frequency-amplitude relations for these clamped and simply supported 
beams with various force intensity (F 0 ) are also plotted in Figures 10 
and 11 for the three-mode solution, respectively. 

The responses of clamped and simply supported immovable beams 
(L/R = 1010) with IDI are shown in Tables 11 and 12. Table 11 shows the 
frequency ratios for the cases of single, two and three-mode 
approaches. It also shows that the higher modes are more important for 
the clamped beam than the simply supported beam. This can be observed 
by looking at one of the amplitude, eg. A/R = + 5.0, the frequency ratio 
for the clamped beam is changed for different numbers of modes used in 
the formulation, but the frequency ratio for the simply supported beam 
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Table 9 Frequency Ratios for Nonlinear Forced Vibration of 

Clamped and Simply Supported Beams (L/R = 1010) without 
Inplane Displacement and Inertia (IDI) under Uniform 
Harmonic Distributed Force 


A/R 

Frequency Ratio, 


1 mode 

2 modes 

3 modes 



Clamped Beam: F Q = 

0.002 N/mm 


- 

1.0 

.4101 

.4105 

(3) a 

.4097 

(3) 



1.3856 

1.3855 

(3) 

1.3856 

(2) 

• 

2.0 

.8592 

.8595 

(4) 

.8573 

(3) 



1.2705 

1.2701 

(3) 

1.2694 

(4) 

± 

3.0 

1.0509 

1.0511 

(5) 

1.0440 

(5) 



1.2994 

1.2987 

(4) 

1.2940 

(4) 

± 

4.0 

1.2189 

1.2189 

(6) 

1.2019 

(7) 



1.3869 

1.3860 

(6) 

1.3725 

(6) 

± 

5.0 

1.3877 

1.3878 

(8) 

1.3554 

(9) 



1.5087 

1.5078 

(8) 

1.4794 

(8) 



Simply Supported Beam: F 0 = 

0.001 N/mm 


- 

1.0 

1.8328 

1.8331 

(3) 

1.8331 

(3) 

± 

2.0 

.8150 

.7937 

(5) 

.7934 

(5) 



1.6840 

1.6771 

(5) 

1.6771 

(5) 

± 

3.0 

1.4013 

1.3560 

(7) 

1.3549 

(7) 



1.8470 

1.8168 

(7) 

1.8162 

(7) 

± 

4.0 

1.8593 

1.7760 

(11) 

1.7725 

(11) 



2.1314 

2.0629 

(11) 

2.0603 

(11) 

± 

5.0 

2.2920 

2.1625 

(17) 

2.1543 

(18) 



2.4742 

2.3581 

(17) 

2.3511 

(18) 


a. Number in brackets denotes the number of iterations to get a 
converged solution 10”^ 


Table 10 Amplitude Ratios for Nonlinear Forced Vibra- 
tion of Clamped and Simply Supported Beams 
(L/R * 1010) without Inplane Displacement 
and Inertia ( ID I ) under Uniform Harmonic 
Distributed Force 


A/R 


Amplitude Ratio 


2 modes 

3 modes 


Aj/A 2 

Aj/A 2 A1/A3 



Clamped 

Beam: F 0 = 0.002 N/mm 


± 

1.0 

- 1054 

- 1056 

2164 

± 

2.0 

- 280 

- 282 

559 

± 

3.0 

- 137 

- 139 

261 

± 

4.0 

- 87 

- 89 

157 

i 

5.0 

- 64 

- 66 

109 

Simply Supported Beam: F 0 = 0.001 N/mm 

- 

1.0 

446 

446 

★ 

± 

2.0 

125 

125 

7652 

± 

3.0 

66 

66 

1972 

± 

4.0 

45 

44 

848 

± 

5.0 

35 

34 

480 

* 

Number 

is larger than 

10 4 . 
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Table 11 Frequency Ratios for Nonlinear Forced Vibration 
of Clamped and Simply Supported Immovable Beams 
(L/R = 1010) with Inplane Displacement and 
Inertia (IDI) under Uniform Harmonic Distributed 
Force 


A/R 


Frequency Ratio, 


1 mode 

2 modes 

3 modes 



Clamped Beam: F 0 = 0 

.002 N/mm. 


± 

1.0 

.3925 

.3929 (3) a 

.3928 

(3) 



1.3805 

1.3804 (3) 

1.3804 

(3) 

± 

2.0 

.8258 

.8260 (4) 

.8258 

(4) 



1.2481 

1.2478 (3) 

1.2478 

(3) 

± 

3.0 

.9888 

.9888 (4) 

.9881 

(4) 



1.2497 

1.2489 (3) 

1.2486 

(4) 

± 

4.0 

1.1227 

1.1220 (5) 

1.1205 

(5) 



1.3033 

1.3017 (4) 

1.3007 

(4) 

± 

5.0 

1.2546 

1.2527 (6) 

1.2497 

(6) 



1.3872 

1.3845 (6) 

1.3821 

(6) 



Simply Supported 

Beam F Q = 0 

.001 N/mm 


- 

1.0 

1.8156 

1.8156 (2) 

1.8156 

(2) 

± 

2.0 

.6436 

.6436 (2) 

.6436 

(2) 



1.6080 

1.6080 (2) 

1.6080 

(2) 

± 

3.0 

1.1837 

1.1837 (2) 

1.1837 

(2) 



1.6879 

1.6879 (2) 

1.6879 

(2) 

± 

4.0 

1.5675 

1.5675 (2) 

1.5675 

(2) 



1.8823 

1.8823 (2) 

1.8823 

(2) 

± 

5.0 

1.9211 

1.9211 (2) 

1.9211 

(2) 



2.1352 

2.1352 (2) 

2.1352 

(2) 

a. 

Number 

in brackets denotes the number 

of iterations 

to 


get a 

converged solution 

10 -5 




Table 12 

Amplitude Ratios 
Clamped Immovable 
Displacement and 
Distributed Force 

for Nonlinear Forced Vibration of 
Beam (L/R = 1010) with Inplane 
Inertia ( ID I ) under Uniform Harmonic 
; F 0 = 0.002 N/mm. 



Amplitude Ratio 


A/R 

2 modes 

3 

modes 


Aj/A£ 

Aj/A2 

Aj/Aj 

• 1.0 

- 1354 

- 1354 

* 

± 2.0 

- 346 

- 346 

2689 

± 3.0 

- 159 

- 159 

1202 

± 4.0 

- 94 

- 94 

682 

± 5.0 

- 64 

- 64 

442 

* Number 

is larger than 10 

4 
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apparently remains the same no matter how many modes are used in the 
formulation. Figure 12 shows the frequency-amplitude relation for the 
clamped immovable beam under uniform harmonic force intensity of F 0 = 0 
(free-vibration case), 0.002 and 0.004 N/mm. It should be noted that 
all curves in this figure shows the hardening type nonlinearity. Figure 

13 shows the comparison of harmonic balance method 17 , experiment 18 and 
finite element method for a clamped immovable beam under uniform 
harmonic distributed force intensity F Q = 0.004170277 N/mm. It clearly 
demonstrates the remarkable agreement between the experiment and the 
finite element (with IDI) solution. 

The responses for the movable cases are shown in Table 13 and 14. 
Table 13 shows the frequency ratios and amplitude ratios for three-mode 
clamped movable beam with IDI under uniform harmonic force. Similarly, 
the results of the simply supported beam are shown in Table 14. Figure 

14 shows the frequency-amplitude relation for three-mode clamped movable 
beam of slenderness ratio L/R = 20. All of the curves in this figure 
shows that the beam eventually exhibits slightly softening type 
nonlinearity. Figure 15 shows the comparison between the immovable case 
and the movable case for a three-mode clamped beam (L/R * 100) with 
IDI under the force intensity F 0 = 20 N/mm. This clearly shows that the 
movable case reduces beam nonlinearity when compared to the immovable 
case. 


5.4 Nonlinear Response to Concentrated 
Harmonic Force 

The application of the finite element method to simulate the case 
of a concentrated force is to let the length of the loaded element 




RMS A/h 



Table 13 Frequency Ratios and Amplitude Ratios for 
Three-Mode Nonlinear Forced Vibration of 
Clamped Movable Beam with Inplane Displace 
ment and Inertia { ID I ) under Uniform Harmonic 
Distributed Force 


A/R 

Frequency 

Amplitude Ratio 

Ratio, 




w NLl /w Ll 

A 1 /A 2 

A l /A 3 


L/R = 100; F 0 

= 20 N/mm 


± 1.0 

.4011 (2) a 

★ 

★ 


1.3560 (2) 

★ 

★ 

± 2.0 

.7617 (2) 

★ 

★ 


1.1911 (2) 

★ 

* 

± 3.0 

.8482 (2) 

* 

★ 


1.1306 (2) 

★ 

★ 

± 4.0 

.8880 (2) 

★ 

★ 


1.0987 (2) 

8271 

★ 

± 5.0 

.9107 (2) 

7704 

★ 


1.0788 (2) 

5488 

★ 


L/R = 20; F 0 = 

5000 N/mm 


± 1.0 

.8137 (3) 

9653 

★ 

' 

1.1538 (2) 

4797 

★ 

• 2.0 

.9061 (3) 

1940 

* 


1.0734 (3) 

1380 

* 

• 3.0 

.9282 (3) 

817 

* 


1.0387 (2) 

651 

- 8020 

± 4.0 

.9319 (3) 

452 

- 5620 


1.0137 (3) 

381 

- 4757 

± 5.0 

.9268 (3) 

290 

- 3655 


.9913 (3) 

252 

- 3202 


a. Number in brackets denotes the number of iterations to 
get a converged solution 10" 5 



Table 14 Frequency Ratios and Amplitude Ratios for Three- 

Mode Nonlinear Forced Vibration of Simply Supported 
Movable Beam with Inplane Displacement and Inertia 
(IDI) under Uniform Harmonic Distributed Force 


A/R 

Frequency 

Ra ti o 

“NLl /u Ll 

Amplitude Ratio 
A1/A2 A1/A3 


L/R = 100; F 0 

= 3 N/mm 


± 1.0 

.6131 (2) 

★ 

* 


1.2744 (2) 

★ 

★ 

± 2.0 

.8293 (2) 

★ 

★ 


1.1453 (2) 

★ 

★ 

± 3.0 

.8897 (2) 

★ 

* 


1.0988 (2) 

* 

★ 

± 4.0 

.9183 (2) 

★ 

* 


1.0747 (2) 

★ 

* 

± 5.0 

.9349 (2) 

★ 

* 


1.0598 (2) 

* 

* 


L/R = 20; F 0 = 1000 

N/mm 


± 1.0 

.8162 (2) 

* 

★ 


1.1537 (2) 

* 

* 

± 2.0 

.9105 (2) 

* 

* 


1.0770 (2) 

9929 

★ 

± 3.0 

.9371 (2) 

5825 

★ 


1.0475 (2) 

4660 

★ 

± 4.0 

.9471 (3) 

3203 

★ 


1.0295 (2) 

2709 

★ 

± 5.0 

.9499 (3) 

2034 

* 


1.0154 (2) 

1779 

★ 


a. Number in brackets denotes the number of iterations to 
get a converged solution 10”^ 
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FREQ 



IMMOVABLE 



FREQUENCY RATIO oj 
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become smaller and smaller. By letting l to be the length of loaded 
element from coordinate x=a to x=b, Eq. (3.55) can be applied to 
evaluate the constant c which provides the harmonic force matrix [h] in 
Eq. (3.59). The effect of the length of the loaded element, A q , is 
studied and shown in Table 15 for a three-mode clamped immovable beam 
(L/R = 1010) with IDI for the total force P of 0.3 N at the middle of 
beam. The simulated distributed force intensity over the loaded element 
is calculated by F Q * P/A for Eq. (3.59). Similarly, the effect 

of i Q is shown in Table 16 for a two-mode simply supported immovable 
beam (L/R = 1010) with IDI. The comparison of a beam under the same 
amount of the total force P a 0.3 N, but different kind of loading, is 
shown in Figure 16. In this figure, a three-mode clamped immovable beam 
(L/R = 1010) with IDI under a concentrated force at the middle 
(A q /l s 2%) is plotted against a similar beam under uniformly distributed 
force over the entire beam (F Q = 0.002 N/mm.). Similarly, the two-mode 
simply supported solutions for the total force P = 0.15N, F 0 = 0.001 
N/mm for uniform distributed force over entire beam case, are plotted in 
Figure 17. It shows that the concentrated force cases are much more 
severe than the uniform distributed force for the cases studied. 



Table 15 Frequency Ratios for Three-Mode Forced Vibration 

of Clamped Immovable Beam (L/R = 1010) with Inplane 
Displacement and Inertia ( ID I ) under Concentrated 



Harmonic Force: 

Total Force P = 0.3 N 



Frequency Ratio, 


A/R 


U q /L) percent 



5 

2 

1 

- 1.0 

1.6425 (3) a 

1.6436 (3) 

1.6437 (3) 

± 2.0 

.5372 (4) 

.5357 (4) 

.5354 (4) 


1.3965 (4) 

1.3971 (4) 

1.3972 (4) 

± 3.0 

.8468 (5) 

.8462 (5) 

.8461 (5) 


1.3485 (4) 

1.3489 (4) 

1.3490 (4) 

± 4.0 

1.0314 (6) 

1.0310 (6) 

1.0310 (6) 


1.3724 (5) 

1.3727 (5) 

1.3728 (5) 

± 5.0 

1.1879 (7) 

1.1876 (7) 

1.1876 (7) 


1.4356 (5) 

1.4359 (5) 

1.4359 (5) 


a. Number in brackets denotes the number of iterations to 
get a converged solution 10" 5 
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Table 16 Frequency Ratios for Two-Mode Nonlinear Forced 

Vibration of Simply Supported Immovable Beam (L/R = 
1010) with Inplane Displacement and Inertia (IDI) under 
Concentrated Harmonic Force: Total Force P = 0.15 N. 


A/R 


Frequency Ratio, w NL1 

/w Ll 


(A 0 /L) percent 


5 

2 

1 

- 1.0 

2.1286 (2) a 

2.1296 (2) 

2.1297 (2) 

- 2.0 

1.7897 (2) 

1.7903 (2) 

1.7904 (2) 

± 3.0 

.9949 (2) 

.9941 (2) 

.9940 (2) 


1.8056 (2) 

1.8060 (2) 

1.8061 (2) 

± 4.0 

1.4658 (2) 

1.4654 (2) 

1.4654 (2) 


1.9625 (2) 

1.9628 (2) 

1.9628 (2) 

± 5.0 

1.8558 (2) 

1.8555 (2) 

1.8555 (2) 


2.1923 (2) 

2.1925 (2) 

2.1925 (2) 


a. Number in brackets denotes the number of iteratives to get 
a converged solution 10" 5 






concentrated force vs. distributed force. 
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Chapter 6 
CONCLUSIONS 

Multiple-mode nonlinear forced vibration of a beam has been 
analyzed by the finite element method. Inplane (longitudinal) 
displacement and inertia (IDI) are considered in the formulation. 
Nonlinear free vibration can be simply treated as a special case of the 
general forced vibration by setting the harmonic force matrix equal to 
zero. The effect of the higher modes is more pronounced for the clamped 
supported beam than the simply supported one. In addition, the beam 
without inplane displacement and inertia yields more effect in the 
higher modes than the one with inplane deformation and inertia. The 
effect of midplane stretching due to large deflection is to increase the 
nonlinearity. However, the effects of inplane displacement and inertia 
are to reduce nonlinearity. A beam with end supports restrained from 
axial movement (immovable case), only the hardening type nonlinearity is 
observed. For beams of large slenderness ratio (L/R > 100) with movable 
end supports, the increase in nonlinearity due to large deflection is 
partially compensated by the reduction in nonlinear due to inplane 
displacement and inertia. This leads to a negligible hardening type 
nonlinearity, therefore, the small deflection linear solution can be 
employed. However, beams of small slenderness ratio (L/R = 20) with 
movable end supports, the softening type nonlinearity is found. The 
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concentrated force case yields a more severe response than the uniformly 
distributed force case. 

The finite element method, in practice, is more suitable to use to 
analyze modern complex structures. Nonlinear theory can be employed to 
obtain more accurate solutions and explain new phenomena. By combining 
the finite element method and nonlinear theory together, the more 
realistic models of structural responses can be resulted easily and less 
time consuming. The nonlinear finite element method which is studied 
herein, may be extended to study many more advanced topic researches, 
for example, the service life of the structure (S-N curve), the study in 
nonlinear random vibration and the effects of sub or super-harmonic 
excitations. This clearly demonstrates an advantage for researchers to 
employ this multiple-mode nonlinear vibration in the future. 
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APPENDIX 


CONVERGENCE CRITERIA 


Three displacement convergence criteria (norms) used by Bergan and 
Clough*^ for multiple-mode nonlinear free and forced vibrations by 
finite element method are employed. These three norms are the maximum 
norm, the modified absolute norm and the modified Euclidean norm. 

The maximum norm is defined as 


e|| M = max 
J 


Av o 


v j,ref 


The modified absolute norm is defined as 






V j,ref 


The modified Euclidean norm is defined as 


(A.l) 


(A. 2) 




V j,ref 


2 

) 


1/2 


(A. 3) 


In these expressions, Av. is the change in displacement component j 

J 

during iterative cycle n, and v^. is the reference displacement which 
is the largest displacement component of the corresponding "type". For 
instance in a nonlinear vibration problem involving deflections w and 
rotations w,^, the reference displacement are the largest deflection and 
the largest rotation of the corresponding type. 


